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D Abstract 
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■4— > Under pressure, the itinerant helimagnet MnSi displays unusual magnetic properties. We have previously discussed a BCC helical 
^ spin crystal as a promising starting point for describing the high pressure phenomenology. This state has topologically nontrivial 
configurations of the magnetization field. Here we note the consequences for magneto-transport that arise generally from such spin 
(— j textures. In particular a skyrmion density induced 'topological' Hall effect, with unusual field dependence, is described. 
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1. Introduction 

Recently, the interplay between electron transport and 
magnetism has been at the center of renewed activity. The 
anomalous hall effect in ferromagnets is a prototypical ex- 
ample of this interplay. Furthermore, transport in the pres- 
ence of more complicated magnetic order, such as spin spi- 
rals, are proving to be extremely rich. An interesting analog 
of the anomalous Hall effect of ferromagnets, the 'topolog- 
ical' Hall effect, is expected [1] where Berry's phase from 
skyrmion spin textures leads to an effective magnetic flux. 
Experimentally, novel physical phenomena have also been 
reported in itinerant helimagnets (with spiral magnetic or- 
der), especially in the B20 structure family which includes 
FeGe and MnSi. The latter was long regarded as an example 
of a well understood magnetic system, until recent experi- 
ments under moderate pressure revealed a quantum phase 
transition into a 'partial order' state, with an unusual signa- 
ture in neutron scattering and non-Fermi liquid resistivity 
signature Ap(T) oc T 15 [2,3]. This non- Fermi liquid resis- 
tivity is seen for almost three decades in temperature from 
around 6K down to a few mK. While the unconventional 
magnetism disappears on further increasing the pressure, 
the non-fermi liquid signatures in transport persists up to 
the highest pressures studied. Similar phenomena are seen 
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in the related compound FeGe, where spin spirals are also 
believed to play a role [4] . 

Recently, the anomalous Hall effect has been measured 
by several groups both at ambient and higher pressures [5] . 
M. Lee et al. have identified a sharp feature in the Hall 
resistivity, which appears to be increased in a restricted 
parameter region of temperature, pressure and field [6]. 
They have put forward the interesting speculation that this 
feature may be related to a 'topological' Hall effect (or Hall 
effect produced by spin chirality) [1]. 

Here, we first review our proposal of BCC helical spin 
crystals as a starting point to describing the 'partial or- 
der' state [7,8]. These and other spin crystals may be un- 
derstood as periodic arrangements of skyrmion configu- 
rations. The physical manifestations of such topologically 
non-trivial spin textures is the main subject of this arti- 
cle. In particular, we discuss the possibility of observing 
a topological Hall effect in such textures, form two simple 
examples. 

Phenomenological theory of anomalous magnetism in 
MnSi at high pressure: MnSi is found to have the fol- 
lowing phase diagram. At ambient pressure, a ferromag- 
netic ordering occurs below 30K. The presence of weak 
Dzyaloshinskii-Moriya interactions which are allowed by 
the non-centrosymmetric B20 cubic structure, leads to a 
spiraling of the magnetization. A single spiral state with 
a pitch of about A « 170A and a small moment 0.3/is 
is produced. Crystalline anisotropy locks the spiral along 
certain fixed crystal directions - a simple symmetry anal- 



ysis reveals that the simplest energy terms consistent 
with the B20 crystalline symmetries are such that they 
favor the spiral along the (111) direction or the (100) di- 
rections. Indeed, neutron scattering at ambient pressure 
reveals Bragg spots along the (111) direction (where (111) 
refers to the [111] and the seven other crystallographically 
equivalent directions). On increasing the pressure beyond 
p c = \A.6kBar, there is a quantum phase transition into 
an unusual state termed the 'partial order' state in Ref. 
[3]. Here, in contrast to ambient pressure, the elastic neu- 
tron scattering reveals enhanced scattering all around the 
wavevector sphere k = 2tt/X. However the intensities are 
peaked around the (110) directions. It may be argued that 
this is very unlikely to be due to single spiral states [7,8]. 

We therefore adopt the next simplest hypothesis, that a 
useful starting point to explain partial order is a spin crys- 
tal, made of a coherent superposition of magnetic spirals 
along the (110) directions. Approaches which are comple- 
mentary to ours Refs. [9,10], have proposed different peri- 
odic spin textures starting from real space rather momen- 
tum space. At longer scales, magnetic order is believed to 
be destroyed by disorder or fluctuation effects. Experimen- 
tally, while NMR sees a frozen magnetic moment in this 
regime [11], yuSR does not [12]. Here, we assume the order 
is static, at least in the presence of a magnetic field. 

Elsewhere, we have studied the energetics of such a state 
relative to the single mode state [7,8]. Spatially varying in- 
teractions between spirals are required to stabilize such a 
state, and is more likely the shorter the spiral wavelength. 
In this respect, although MnSi has a relatively long spi- 
ral wavevector, it is the shortest compared to others in 
the same family (e.g. FeGe, FeCoSi etc.). Here we do not 
revisit the energetics, but rather adopt a phenomenolog- 
ical approach and deduce the consequences of having a 
(110) spin crystal. Note, the most general such spin crys- 
tal can be specified by the amplitude and phase of the spi- 
rals along each direction. This requires specifying six com- 
plex numbers tjji , ip2 ■ ■ ■ i>§ ■ From these, and a specification 
of the basis used along each spiral direction, one can re- 
construct M(r). If we choose the six wavevectors k, to be 
along {(1, 1, 0), (1, 1, 0), (0, 1, 1), (0, 1, 1), (1,0, 1), (1, 0, !)}, 
and the basis vectors to be unit vectors e'j oc [(0, 0, 1) x kj], 
and e'j oc [e" x kf], we can write: 

M(r) = i]TV^+^ k - r + ^- (1) 

These spin structures have the periodicity of a body cen- 
tered cubic (bec) lattice, hence, bec spin crystals. 

We can now develop a Landau theory in terms of these 
complex fields. Cubic symmetry restricts the Landau free 
energy to the following form: 

F = E V H' M 2 Wi' I 2 + A Rc ( T * +Ty + T z ), (2) 
i? 

where T x = ViV^V'sV'e, T y = tptip^ipi and T z = 
— V'3V'4V'5'06- The bec state is realized if the parameters Vij 



and A are such that the energy minimum has equal ampli- 
tudes for all six modes. As for the phases, it is clear, that 
they only depend on the sign of the parameter A. For A > 
(< 0), the phases arrange in such a way that T x ,T y ,T z 
are negative (positive). The remaining three phase degrees 
of freedom correspond to global translations of the struc- 
ture. Thus symmetry considerations alone lead to just 
two possible bec states, called beel and bcc2 respectively. 
Both states are doubly degenerate due to time reversal 
M — > -M. The state beel (A > 0) is more likely to be 
realized than bcc2, because the magnitude of magnetic 
moments |M(r)| has a much narrower distribution in beel 
than in bcc2 [7]. The resulting real-space pattern M(r) of 
beel is best understood as a periodic network of skyrmion 
lines and anti- vortex node-lines [7,8]. These structures are 
stabilized by symmetry - the magnetic point group sym- 
metry of the structure is 0(T). This group is generated 
by 7r/2 rotations around the x-, y- or z-axis, followed by 
time- reversal [7,8]. We will see that the most direct man- 
ifestation of these structures is in the magneto-transport, 
leading to a 'topological' Hall effect. 

2. Skyrmion density 

Non-trivial magnetic structures in metals can act as an 
effective orbital magnetic field on conduction electrons and 
thus lead to an anomalous Hall effect which is indepen- 
dent of spin-orbit coupling [1]. Essentially, when conduc- 
tion electrons are forced to follow the local spin direction 
M(r), they acquire a Berry's phase. This may be viewed 
as an effective field, which is proportional to the skyrmion- 
density <&(r), defined as 

= -J-e^ n • (dph x <9 7 n) , (3) 

87T 

where n = M/|M|. The flux of $ through any two- 
dimensional rectangle with periodic boundary conditions 
is a topological winding number, i.e. an integer. Note, the 
terminology 'spin chirality' is somewhat misleading since 
$ transforms as an axial vector, like the magnetic field. 
It is even under inversion, hence not a chiral object - i.e. 
broken inversion symmetry is not required to have $ =/= 0. 

It is remarkable that so many different theoretical pro- 
posals for novel magnetic structures in 3D magnets without 
inversion symmetry independently lead to arrangements 
of skyrmion lines (meron lines, double-twist cylinders) 
[7,8,9,10]. Therefore, such systems appear to be natural 
candidates for observing the effects of <E>. 

By inspection of Eq. (3), it becomes immediately clear 
that a simple spin density wave, which involves a single pair 
of wavevectors ±k, always has a vanishing skyrmion den- 
sity, since in this case M depends only on a single space co- 
ordinate. This is also true if M acquires a uniform compo- 
nent, i.e. in the conical phase under a magnetic field. Multi- 
q spin crystals are therefore the simplest slowly varying 
magnetic structures with non- vanishing skyrmion density. 




Fig. 1. The dotted line shows the skyrmion number J dxdy<& z of 
one unit cell for the helical 2D square spin crystal with a uniform 
component M(r) = (sin y, sinrr, m + cos y — cos x) as a function of 
m. The full line shows the same quantity evaluated with a cutoff for 
M 2 < 1/3, as explained in the text. 

2.1. Square spin crystal 

As the simplest example for helical spin crystal, we first 
consider the magnetic texture obtained by superimposing 
two orthogonal spirals of the same chirality. For simplicity, 
we choose units such that each spiral has unit amplitude 
and a wave-length equal to 2ir. In appropriate coordinates, 
the square spin crystal is then given by 

M(r) = (sin y, sin x, cosy — cos x) (4) 

Note, it is a 2D like structure given that there is no vari- 
ation along the z direction. This is a periodic 2D arrange- 
ment of skyrmion-like textures pointing up and down in 
a checkerboard pattern. The skyrmion density <I> Z follows 
this pattern, in fact, $ z = M z /(47r|M| 3 ) in this case. Thus 
for symmetry reasons, the total skyrmion number over the 
unit cell (fdxdy$ z ) vanishes, since positive and negative 
contributions cancel exactly. Now, a magnetic field is ap- 
plied in the z-direction. We assume, in the interest of clar- 
ity, that this only adds a uniform magnetization (0, 0, m) 
to M, but does not otherwise affect the spin arrangement. 
The skyrmion number per unit cell is an integer, and there- 
fore can not change continuously with m. The behavior is 
shown in Fig. 1. The function is discontinuous at m = 0, ±2 
since it is at precisely these values of m that the spin tex- 
ture develops nodes. These values can be easily obtained 
from the m = magnetization pattern of Eq. 4 by noting 
that the magnetization that is purely along the z axis takes 
on the values M = 0, ±2z. In the parameter regime < 
\m\ < 2, a non-trivial, quantized spin chirality is obtained, 
$ z = — sgn(m). 

At least in the limit of strong coupling between conduc- 
tion electrons and magnetic background, it is natural that 
the spin direction n rather than the amplitude enters Eq. 
(3). However, in contrast to Refs. [1], we are concerned with 
a smoothly varying magnetization M(r), which has a spa- 
tially varying amplitude and even nodes. In the vicinity of 
nodes, electrons can not be infinitely sensitive to directional 
changes of the magnetic background and the theory needs 
to be extended. To include this effect in a simple way, let 
us assume that the electrons couple only then to the back- 
ground spin texture, if the latter has an amplitude bigger 
than some cutoff. I.e. we consider $ c = 6 (M 2 - M§) $, 



where is the Heaviside step function. As a result of 
this modification, the discontinuous jumps of the skyrmion 
number gets smoothed out, as shown in Fig. 1. 

2.2. bccl spin crystal 

We now consider the spin crystal bccl. An explicit ex- 
pression for M(r) is given in Ref. [8]. Again, the average 
of the skyrmion density over a unit cell is zero for symme- 
try reasons [the point group symmetry of bccl is 0(T) and 
therefore does not tolerate a nonzero axial vector like ($)]. 
Hence, there are local skyrmion densities and conduction 
electrons will feel their effective fields, but the average ef- 
fective field is zero. No anomalous Hall effect can therefore 
occur without an applied field. 

Applying an external magnetic field lowers the symme- 
try and allows for ($) ^ 0. An external field has two effects 
on the bccl structure. First, it will induce a uniform mag- 
netization m pointing along the field and reduce the am- 
plitude of the helical components. Second it will affect the 
relative amplitude and phases of the six interfering helical 
modes, thus transforming the spin texture. Eventually, the 
magnetic field will induce a transition to a different mag- 
netic state, for example the conical single-spiral state or 
the spin-polarized (ferromagnetic) state. Here, we will only 
discuss the effect of a uniform component m added to the 
otherwise unperturbed bccl structure. 

Before we go into the details, we make a short remark 
about time reversal symmetry (T). The bccl state breaks 
T in such a way that it may not be restored by a subsequent 
translation [7,8,13]. As a consequence, there are two degen- 
erate bccl states related by time reversal: M(r) — > — M(r). 
Let us label them by a parameter S = ±1. Time reversal 
symmetry requires that the integrated skyrmion density of 
the state S plus a uniform component m satisfies 

<*>U,m = -<*>l-S,-m (5) 

In the following, we choose the lattice constant of the 
cubic (i.e. not elementary) unit cell as the unit of length and 
the amplitude of the elementary spiral modes which build 
up bccl as the unit of magnetization. The net effective field 
is determined by the integral of <E> over the unit cell. It may 
be trivially written as 

($ z ) - Jdzn z (z) (6) 

where n z = J dxdy $ is an integer- valued function of z. It 
can only change its value at those positions of z, which 
correspond to nodes of M. For m = 0, the bccl structure 
has node lines, and n z (z) could potentially change every- 
where. Fortunately, we know from symmetry that in this 
case ($) = 0. For m ^ 0, there are only discrete point 
nodes, and Eq. (6) turns out to be very useful to evaluate 
($ z ). The other components, ($ x ' v ) may be obtained in a 
similar way. 
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Fig. 2. Spin Chirality for bccl structure with a uniform magnetization 
(0, 0, m) superimposed (a) Shows (& z ) as function of m. (b) Same 
with a cutoff for M 2 < 1. 

As an example, we consider the case m = (0, 0, m), i.e. 
a uniform component along z is added to the bccl tex- 
ture (the case with the field along (111) direction is found 
to be qualitatively similar). For this geometry, ($) is an 
odd function of m and parallel to 2, because of the 0(T) 
point group symmetry. Hence, by Eq. (5), the result is in- 
dependent of the time-reversal label S. The nodes of M are 
particularly easy to analyze as a function of m. As a re- 
sult, (& z ) may be calculated analytically. We find (<£> 2 ) = 
-sgn(m) (2/7r) arccos(|m|/2/V2) for < \m\ < 2\/2 and 
($ 2 ) = elsewhere. The behavior is shown in Fig. 2. For 
comparison, the maximum spin density |M| at m = is 
2\/3. There are two phase transitions as a function of the 
uniform component. One is of course at m = 0, where 
($ 2 ) changes sign. However, the nontrivial transition oc- 
curs at the high-field ends, where ($ z ) vanishes continu- 
ously but with an infinite slope. Note, at this point the mag- 
netic structure is still non-trivial, but the 'spin-chirality' 
has been ironed out by the large uniform magnetization. 
Like in the former example, introducing a cutoff magnitude 
Mo, to model the decoupling of electrons from the low mag- 
netization parts of the textures, smoothes these transitions 
to some degree. For the BCC1 spin crystal, with a spiral 
wavelength of 170A, the 'spin-chirality' induced effective 
field is about 1.7Tesla. 

3. Summary 

We have investigated the possibility of a topological 
Hall effect in non-trivial magnetic structures composed of 
skyrmion configurations under the influence of an external 
magnetic field. In two simple examples, we have analyzed 
the integrated skyrmion density ($) , which has been shown 
to act on conduction electrons as an effective internal field. 

In both cases it is clear that an induced magnetization 
plays two opposite roles. On the one hand, it is required by 
symmetry to induce a non-zero spin chirality, but on the 
other hand, increasing it beyond a threshold value leads to a 



vanishing ($) . This will ultimately lead to a non-monotonic 
effect of magnetization on the Hall conductivity, which is 
influenced by (<&) . This, we believe is a hallmark of the spin 
chirality contribution to the Hall effect. 

In the case of a two-dimensional spin structures, which 
may be realized for example in thin films, magnetic surfaces 
or interfaces, There is the possibility of a plateau region, 
where the effect is independent of the external field, due to 
quantization of the skyrmion number. 

In both above cases, the effective field created by spin 
configurations turns out to be antiparallel to the uniform 
magnetization. This predicts the topological Hall effect to 
be of opposite sign to the normal Hall effect. 

Questions left for future investigation include whether 
the 'topological' Hall effect discussed here is relevant to 
the experiments of [6], which are still on the low pressure 
side, and how electron transport is affected by a disordered, 
and possibly fluctuating, arrangement of non-trivial spin 
textures. 
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